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1  Introduction 

The  purpose  of  this  grant  was  to  provide  support  to  allow  a  student  working  under  the 
preceding  grant,  F49620-94-1-0017  to  transition  to  new  funding  and  to  finish  his  graduate 
degree.  Both  objectives  were  met.  The  research  program  involved  two  projects.  The  first 
focused  on  the  development  of  extreme  ultraviolet  (XUV)  sources,  and  their  application 
to  the  development  of  high  resolution,  microscopic  imaging  techniques.  That  project  came 
to  close  shortly  after  the  beginning  of  this  grant.  The  second  project  investigated  unique 
approaches  to  building  multiuser  optical  fiber  networks. 


2  Personnel  Supported  This  Period 

Senior  Personnel 

James  F.  Young,  Principle  Investigator  and  Professor. 


Students 

Tasshi  Dennis,  graduate  student;  Ph.D.  completed  September  1998. 
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Honors 


Tasshi  Dennis  was  awarded  a  National  Research  Council  two  year  Research  Associate  Fel¬ 
lowship  at  NIST,  Optoelectronics  Division,  Fiber  and  Discrete  Components  Group,  Boulder, 
CO. 

3  Summary  of  Research 

3.1  Short  Wavelength  Imaging 

This  grant  provided  supplemental  and  continuing  student  support  for  the  preceding  grant, 
F49620-94-1-0017,  which  ended  14  October  1997.  The  short  wavelength  imaging  project 
was  terminated  at  that  time  because  of  technical  and  funding  reasons.  Consequently,  all 
effort  on  this  grant  in  that  area  has  been  joint  with  F49620-94-1-0017,  and  the  reader  is 
referred  to  its  Final  Technical  Report  for  details.  The  final  publication  in  the  area  is  listed 
below,  and  attached  to  this  report. 

3.2  Fiber  Optic  Networks 

The  bulk  of  the  effort  supported  by  this  grant  focused  on  multi-user,  multi-wavelength 
optical  fiber  communications  systems.  Wavelength  division  multiplexing  (WDM)  is  a  very 
successful  technique  for  increasing  the  capacity  of  point-to-point  fiberoptic  communications 
links.  We  have  been  investigating  more  sophisticated  encoding  techniques  that  promise 
even  further  increases  in  functionality,  particularly  for  multiuser  networks.  For  example, 
radio  frequency  (RF)  code  division  multiple  access  (CDMA)  has  been  very  successful  at 
supporting  a  large  pool  of  asynchronous  network  subscribers  sharing  a  common  channel, 
increasing  noise  tolerance,  and  providing  user  security  and  flexibility.  This  success  has 
been  primarily  due  to  the  development  and  use  of  bipolar,  {-1-1,  —1},  code  sequences  and 
bipolar  correlation  detection.  While  both  time  and  spectral  domain  optical  CDMA  have 
been  studied,  the  results  to  date  have  been  poor  because  direct  detection  optical  channels 
are  inherently  unipolar,  {-f-1,0},  and  cannot  directly  implement  the  RF  bipolar  coding  and 
correlation  detection. 

We  have  invented  a  technique  to  achieve  true  bipolar  encoding  and  decoding  for  optical 
spectral  CDMA  [“All-Optical  CDMA  with  Bipolar  Codes,”  L.  Nguyen,  B.  Aazhang,  and 
J.  F.  Young,  Electronics  Letters  31,  469-470  (1995);  U.S.  Patent  No.  5,760,941],  and  per¬ 
formed  proof-of-principle  measurements.  We  have  examined  the  practical  implementation  of 
passive,  all-optical  encoders  and  decoders  incorporating  these  principles,  using  a  fiber-based, 
three  user  testbed  capable  of  OC-24  (1.25  Gb/s)  operation,  and  developed  a  numerical  sim¬ 
ulation  of  the  encoder  and  decoder  correlation  process  to  study  the  various  optical  effects 
that  influence  the  detection  performance.  Recently,  we  have  completed  bit  error  rate  mea¬ 
surements  under  a  variety  of  conditions  using  our  testbed.  Details  of  this  work  can  be  found 
in  the  papers  listed  in  the  next  section. 
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Our  work  has  shown  that  passive,  all-optical  correlators  are  attractive  for  high-speed 
network  applications  where  asynchronous  operation  and  independence  of  data  rate  are  de¬ 
sired.  They  are  able  to  process  the  extremely  high  combined  data  rate  on  a  common  network 
channel  while  the  detection  electronics  operate  at  only  individual  user  rates.  However,  opti¬ 
cal  distortions  cause  degradation  of  correlator  performance  and  have  a  direct  impact  on  the 
capacity  of  a  network  to  support  simultaneously  active  users.  Our  numerical  simulations, 
analytical  results,  and  experimental  measurements  have  identified  wavelength  scale  distor¬ 
tion  as  a  critical  optical  parameter  in  the  correlation  process  and  illustrated  the  importance 
of  well-matched  physical  characteristics  when  using  multiple  gratings  or  other  dispersive  el¬ 
ements.  Based  on  our  results,  we  have  received  a  grant  to  continue  the  work  from  the  Texas 
Advanced  Technology  Program. 

The  usefulness  of  our  technique  relative  to  WDM  will  ultimately  depend  on  a  tradeoff 
between  the  added  complexity  and  the  added  functionality.  Encoder  and  correlator  designs 
that  minimize  complexity  and  optical  distortions  are  thus  very  important.  We  have  begun  a 
collaboration  with  Dr.  Ray  Chen  at  the  University  of  Maryland,  Baltimore  County,  to  design 
and  fabricate  planar  waveguide  integrated  optical  circuits  that  will  implement  our  encoder 
and  correlator  concepts  on  a  single  photonic  chip.  These  devices  will  improve  performance 
while  minimizing  complexity  and  cost. 


4  Publications  and  Presentations  of  Supported  Work 

“Measurements  of  BER  Performance  for  Bipolar  Encoding  of  an  SFS,”  T.  Dennis  and  J.  F. 
Young,  (In  preparation). 

“Optical  Implementation  of  Bipolar  Codes,”  T.  Dennis  and  J.  F.  Young,  (In  Review;  draft 
attached) . 

“Experimental  Demonstration  of  All-Optical  CDMA  using  Bipolar  Codes,”  Tasshi  Dennis, 
Ph.D.  Dissertation,  Electrical  and  Computer  Engineering  Department,  Rice  University, 
Houston,  Texas  (September  1998). 

“Experimental  Demonstration  of  Bipolar  Codes  for  Optical  Spectral  Amplitude  CDMA  Com¬ 
munication,”  L.  Nguyen,  T.  Dennis,  B.  Aazhang,  and  J.  F.  Young,  J.  Lightwave  Tech¬ 
nology  15,  1647-1653  (1997). 

“Laser-Produced  Plasmas  as  Short  Wavelength,  Incoherent  Optical  Sources,”  J.  F.  Young, 
in  Atomic,  Molecular,  and  Optical  Physics:  Electromagnetic  Radiation,  F.  B.  Dunning 
and  R.  G.  Hulet,  eds.,  ,  Chapter  1  (Academic  Press,  1997). 

“Imaging  Characteristics  of  Poly(methyl  methacrylate)  at  Vacuum  Ultraviolet  Wavelengths,” 
I.  E.  Ferincz,  C.  Toth,  and  J.  F.  Young,  J.  Vac.  Sci.  Technol.  B  15,  828-832  (1997) 
(Attached). 
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“All-Optical  Encoders  and  Decoders  for  Spectral  CDMA  using  Bipolar  Codes,”  T.  Dennis 
and  J.  F.  Young;  in  All-Optical  Networking:  Architecture,  Control,  and  Management 
Issues,  SPIE  (in  press). 

“Wavelength  Encoding  for  Optical  Networks,”  T.  Dennis  and  J.  F.  Young;  Int.  Conf.  on 
Applications  of  Photonic  Technology,  Ottawa,  Canada,  July  1998  (Invited). 

“Wavelength  Encoding  to  Reduce  Four-Wave  Mixing  Crosstalk  in  Multi- Wavelength  Chan¬ 
nels,”  Y.  Guo,  B.  Aazhang,  and  J.  F.  Young;  IEEE  Lasers  and  Electro-Optics  Society 
Annual  Meeting,  San  Francisco,  CA,  November  1997. 

“Demonstration  of  All-Optical  CDMA  with  Bipolar  Codes,”  T.  Dennis,  B.  Aazhang,  and 
J.  F.  Young;  IEEE  Lasers  and  Electro-Optics  Society  Annual  Meeting,  San  Francisco, 
CA,  November  1997. 

5  Relationship  to  Other  Programs  and  Support 

The  imaging  applications  portion  of  this  work  was  jointly  supported  by  the  National  Science 

Foundation;  the  optical  communications  work  was  jointly  supported  by  the  Texas  Advanced 

Technology  Program. 


File:  c :  / j f y/graiit3/AF.56/ af 96-final .  tex. 
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Optical  Implementation  of  Bipolar  Codes  * 

Tasshi  Dennis  and  James  F.  Young  {Fellow,  IEEE) 

Department  of  Electrical  and  Computer  Engineering, 
Rice  University,  MS-366,  PO  Box  1892,  Houston,  TX  77251 
Phone:  713-527-4721;  FAX:  713-524-5237;  young@rice.edu 


Abstract 

The  physical  implementation  of  a  bipolar  encoding  scheme  suitable  for  fiber  op¬ 
tic  networks  is  reported  with  both  experimental  and  theoretical  results.  The  power 
spectrum  of  an  erbium-doped  superfiuorescent  fiber  source  is  encoded,  the  bipolar 
correlations  of  the  codes  are  verified  and  rejection  of  multiple-access  interference  is 
demonstrated  in  a  fiber-based  testbed.  Simulations  of  the  correlation  process  identify 
key  optical  parameters  and  physical  characteristics  important  to  the  design  of  future 
systems.  A  modification  to  the  experimental  testbed,  made  according  to  the  theoretical 
analysis,  results  in  improved  correlation  performance. 


1  Introduction 

Wavelength  division  multiplexing  (WDM)  is  a  very  successful  technique  for  increasing  the 
capacity  of  fiberoptic  communications  links.  WDM  is  the  simplest  form  of  wavelength  encod¬ 
ing,  and  we  have  been  investigating  more  sophisticated  encoding  techniques  that  promise 
even  further  increases  in  functionality,  particularly  for  multiuser  networks.  For  example, 
radio  frequency  (RF)  code  division  multiple  access  (CDMA)  has  been  very  successful  at 
supporting  a  large  pool  of  asynchronous  network  subscribers  sharing  a  common  channel, 
increasing  noise  tolerance,  and  providing  user  security  and  fiexibility.  This  success  has  been 
primarily  due  to  the  development  and  use  of  bipolar,  {+1,  —1},  code  sequences  and  bipo¬ 
lar  correlation  detection.  While  both  time  and  spectral  domain  optical  CDMA  have  been 
studied  [1,  2,  3],  the  results  to  date  have  been  poor  because  direct  detection  optical  channels 
are  inherently  unipolar,  {-t-1, 0},  and  cannot  directly  implement  the  RF  bipolar  coding  and 
correlation  detection. 

We  have  previously  reported  a  technique  to  achieve  true  bipolar  encoding  and  decoding 
for  optical  spectral  CDMA  [4],  and  presented  proof-of-principle  measurements  [5]  .  In  this 

*This  work  was  jointly  supported  by  the  Texas  Advanced  Technology  Program  and  the  U.S.  Air  Force 
Office  of  Scientific  Research. 
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paper,  we  examine  the  practical  implementation  of  passive,  all-optical  encoders  and  decoders 
incorporating  these  principles.  First,  we  present  experimental  results  using  a  fiber-based, 
three  user  testbed  capable  of  OC-24  (1.25  Gb/s)  operation.  Second,  we  develop  a  numerical 
simulation  of  the  encoder  and  decoder  correlation  process  to  study  the  various  optical  effects 
that  influence  the  detection  performance.  Finally,  we  make  modifications  to  our  decoder 
based  on  our  theoretical  results  and  show  an  improvement  in  correlation  performance. 


2  Encoder  and  Decoder  Operation 

Our  basic  scheme  has  been  reported  previously  [4,  5],  and  is  based  on  representing  a  bipolar 
code  as  two  unipolar  codes.  For  each  length-AT  bipolar  code  X,  we  form  a  unipolar  sequence 
U  by  replacing  each  —1  with  0.  The  unipolar  super-code  J  of  length  2N  is  then  formed  as  the 
concatenation  of  U  and  its  complement,  U  ®U.  At  the  decoder,  this  technique  reduces  the 
bipolar  correlation  to  two  unipolax  correlations  that  can  be  performed  using  a  passive  optical 
system  and  intensity  detection.  Our  experimental  apparatus  performs  spectral  encoding 
and  decoding  of  a  superfluorescent  fiber  source  (SFS)  using  a  bulk  grating  and  amplitude 
mask  as  shown  in  Fig.  1.  Shown  in  Fig.  2  is  the  fiber-based  testbed  constructed  for  this 


Figure  1:  The  geometry  for  the  grating  based  encoders.  The  spectral  patterns  of  the  two 
outputs,  J  and  J  are  complements  of  one  another,  and  are  used  to  transmit  a  One  and  Zero 
respectively.  The  decoder  has  exactly  the  same  geometry,  but  the  SFS  input  is  replaced  by 
the  signal  from  the  network  channel,  and  the  outputs  are  focused  onto  a  balanced  receiver. 

series  of  measurements.  We  used  a  SFS  consisting  of  40  m  of  highly  doped  erbium  fiber 
pumped  by  a  150  mW,  980  nm  pigtailed  diode.  About  22  mW  of  randomly  polarized  SFS 
power  was  divided  with  a  fiber  polarization  splitter,  forming  broadband  light  sources  for 
two  encoders.  We  encoded  about  9nm  of  the  bandwidth  centered  around  the  1530  nm  peak 
using  Walsh  codes.  These  codes,  which  have  zero  crosscorrelations  in  theory,  were  used  to 
simplify  comparison  between  theory  and  experimental  results.  Three  super-codes  of  length 
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Code  b{0,1} 


Receiver  Code 

Figure  2:  Experimental  setup  of  the  spectral  amplitude  CDMA  testbed;  showing  superfluo- 
rescent  fiber  source  SFS,  isolator/wavelength  division  multiplexer  IWDM,  fiber  splitters  and 
combiners  BS  and  PBS,  polarization  controllers  PC,  and  reflection  and  transmission  outputs 
R  and  T. 

32,  which  we  shall  denote  as  J^,  Jb,  Jc-,  and  their  complements,  were  etched  on  aluminum- 
coated  glass  substrates  for  each  encoder  and  decoder  in  the  testbed.  The  smallest  etched 
feature  was  200  /xm  wide,  corresponding  to  a  spectral  slice  of  approximately  0.18  nm  of  optical 
bandwidth.  To  simplify  mask  fabrication,  guardbands  between  adjacent  code  elements  were 
not  used.  The  two  outputs  of  each  encoder,  formed  according  to  the  selected  bipolar  code, 
are  ON/OFF  modulated  with  2x1  interferometric  1  GHz  optical  switches  according  to  the 
information  symbol  h.  The  modulated  outputs  are  combined  to  form  the  network  channel 
and  sent  to  the  decoder.  The  mask  of  the  decoder  performs  the  bipolar  correlation  according 
to  the  code  selected  for  receiving,  and  the  output  is  detected  with  a  balanced  receiver  and 
monitored  with  a  spectrum  analyzer. 


3  Correlation  Measurements 

The  correlations  between  the  codes  were  measured  with  a  single  encoder  and  decoder.  With 
the  decoder  adjusted  to  receive  a  particular  code,  the  encoder  mask  was  positioned  in  suc¬ 
cessive  measurements  to  transmit  any  one  of  the  available  code  patterns.  This  procedure 
was  repeated  for  each  of  the  codes  at  the  decoder,  while  the  optical  switch  was  driven  by  a 
sinusoid.  As  shown  in  Fig.  3,  for  the  same  code  ( J^,  Jb,  Jc)  at  encoder  and  decoder  a  large 
autocorrelation  power  was  detected,  one  each  for  positive  and  negative  receiver  currents.  For 
mismatched  codes,  the  crosscorrelation  signals  were  on  average  23.2  dB  below  the  autocor¬ 
relations,  showing  that  a  binary  information  symbol  can  be  recovered  by  a  zero  threshold 
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Figure  3:  Correlations  of  the  codes  measured  on  the  spectrum  analyzer.  For  mismatched 
codes,  the  crosscorrelation  signals  are  an  average  23.2  dB  below  the  autocorrelations 

operation.  The  non-zero  crosscorrelations  and  the  variation  of  the  autocorrelations  are  be¬ 
lieved  to  be  due  to  a  combination  of  mask  fabrication  tolerances  and  the  optical  properties 
of  the  system,  which  are  discussed  in  the  next  section. 

The  multiple  access  interference  (MAI)  rejection  was  evaluated  using  two  encoders  and 
one  decoder.  While  one  encoder  and  decoder  were  adjusted  to  matched  codes,  the  second 
encoder  sent  the  other  codes  as  interference  in  successive  measurements.  Figure  4  summarizes 
the  MAI  measurements,  showing  the  change  in  signals  measured  at  the  decoder  when  the 
interference  signal  was  present.  To  confirm  that  the  observed  rejection  was  not  simply 
frequency  multiplexing,  MAI  measurements  were  performed  with  the  encoders  driven  at  the 
same  and  at  difllerent  modulation  frequencies.  In  both  cases  the  conclusion  was  the  same: 
the  degradation  of  the  desired  signal  was  less  than  1.7  dB,  and  the  decoder  could  reject  MAI. 


4  Simulation  Results 

Despite  the  ideal  orthogonality  of  Walsh  codes,  our  experimental  results  do  not  show  zero 
crosscorrelations  or  complete  interference  rejection.  We  believe  this  nonideal  behavior  is 
largely  the  result  of  optical  effects  occurring  in  the  encoder  and  decoder.  In  previous  work 
[5],  we  reported  the  measurement  of  wavelength  scale  distortion  and  limited  resolution.  Each 
of  these  effects  is  apparent  from  the  measured  spectra  of  Fig.  5,  which  shows  the  1530  nm 
spectral  peak  of  the  erbium-doped  SFS,  an  amplitude  mask  pattern,  and  the  resulting  en¬ 
coded  spectrum  seen  at  the  decoder  mask.  The  encoded  pattern  is  about  5%  wider  (500  nm. 
over  a  10  mm  spatial  field)  than  the  encoding  mask.  The  rounded  chip  shapes  indicate  a 
full-width-at-half-maximum  resolution  of  175  /zm.  Accurate  correlations  also  require  that 
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Figure  4:  The  change  in  signals  measured  at  the  decoder  when  an  interference  signal  is 
present. 

each  code  chip  have  equal  power,  and  the  bar-widths  of  the  masks  were  intentionally  varied 
to  compensate  for  the  nonuniform  spectral  power  density  of  the  SFS.  This  compensation  is 
an  approximation  and  the  inaccuracies  degrade  correlation  performance. 

To  better  understand  the  influence  that  all  these  optical  parameters  have  on  correlation 
performance,  a  numerical  simulation  was  performed  using  the  calculation  tool  Matlab.  Using 
a  signal  processing  approach,  the  optical  processes  were  implemented  independently,  from 
the  formation  of  the  super-codes,  to  the  detection  of  the  correlation  signals  with  a  balanced 
receiver.  The  result  of  simulating  the  encoding  of  supercode  Jc  with  the  measured  values 
of  the  optical  parameters  is  shown  in  Fig.  6,  which  compares  favorably  with  the  measured 
spectrum  of  Fig.  5.  Of  the  optical  parameters  considered,  it  would  be  useful  to  know  which 
has  the  greatest  impact  on  the  correlation  process,  so  that  future  systems,  either  ones  in¬ 
corporating  new  technology  or  modifled  versions  of  the  existing  one,  can  yield  the  greatest 
improvement  in  correlation  performance.  Shown  in  Fig.  7  (a)  is  the  simulated  autocorrelation 
signal  as  a  function  of  the  wavelength  scale  distortion  for  code  Jb-  At  zero  distortion,  the 
autocorrelation  is  reduced  to  58%  of  its  optimal  value  by  reduced  chip  contrast  due  to  limited 
resolution.  For  an  increasing  magnitude  of  distortion,  the  autocorrelation  signal  is  reduced 
further,  with  a  slightly  larger  reduction  for  negative  distortion  because  the  total  encoded 
power  is  also  being  reduced,  while  for  positive  distortion  the  power  is  being  held  constant  by 
the  spatial  width  of  the  decoder  mask.  Plot  (b)  shows  the  dependency  of  the  crosscorrelation 
signal  between  decoder  mask  Jb  and  encoded  pattern  Ja,  which  predictably  gets  worse  as 
the  magnitude  of  the  distortion  increases.  Plot  (c)  shows  the  ratio  of  the  of  the  crosscorre¬ 
lation  curve  (b)  to  the  autocorrelation  curve  (a),  which  gives  a  signal-to-interference  ratio 
(SIR),  expressed  as  a  percentage.  At  zero  distortion,  the  SIR  is  4.9%  or  — 26.2dB,  due  to 
the  combined  effects  of  limited  resolution  and  imperfect  compensation  for  the  nonuniform 
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Figure  5:  Erbium  SFS  spectrum,  encoding  mask,  and  an  encoded  spectrum  as  seen  at  the 
decoder  mask.  The  pattern  is  super-code  Jc  —  (1001100110011001  ©  0110011001100110) 
used  in  the  correlation  measurements. 


Figure  6:  Simulation  of  the  encoding  process  based  on  experimentally  measured  optical 
parameters,  showing  an  encoding  mask  and  a  processed  spectrum.  The  results  resemble  the 
measurements  of  Fig.  5. 


6 


(a)  (b) 


-500  -250  0  250  500  -500  -250  0  250  500 


-500  -250  0  250  500  -500  -250  0  250  500 

Distortion  (microns)  Distortion  (microns) 


Figure  7:  Simulated  correlation  performance  as  a  function  of  the  wavelength  scale  distortion. 
Plots  (a),  (b),  and  (c)  were  calculated  for  a  decoder  mask  of  Jb  with  spectra  of  Ja  and  Jb, 
and  plot  (d)  was  for  a  decoder  mask  of  Jb  and  spectra  of  Jb  and  Jc- 
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spectral  shape.  At  the  experimental  value  of  +500  of  distortion,  the  simulated  SIR  rises 
substantially  to  15.6%  or  -16.1  dB,  clearly  demonstrating  the  devastating  impact  of  this 
parameter.  For  comparison,  the  measured  SIR  for  decoder  mask  Jb  and  encoder  mask  Ja 
from  Fig.  3  was  6.5%  or  —23.7 dB.  The  measurements  are  better  than  the  simulation  results 
because  of  the  ambiguity  in  horizontally  positioning  the  decoder  mask  relative  to  the  mask 
at  the  encoder.  For  the  autocorrelation  measurements,  the  correct  horizontal  position  was 
easily  found  by  maximizing  the  received  signal.  For  the  crosscorrelations,  the  ideal  position 
should  be  where  the  received  signal  is  minimum.  However,  horizontal  positioning  during 
a  crosscorrelation  measurement  can  also  compensate  for  the  signal  imbalance  that  occurs 
from  optical  effects  such  as  distortion  and  limited  resolution,  or  even  mask  inaccuracies  and 
defects,  such  as  pin-holes.  A  number  of  alignment  features  on  the  masks  help  to  minimize 
this  problem,  but  the  vertical  translation  of  the  masks  required  to  change  the  codes  adds 
variation  to  an  already  sensitive  alignment. 

Simulations  with  other  codes  at  encoder  and  decoder  revealed  that  the  impact  of  the 
wavelength  scale  distortion  is  also  code-dependent,  as  shown  in  Fig.  7  (d).  The  SIR  curve 
shown  is  for  a  decoder  mask  of  Jb  and  spectra  of  Jb  and  Jc-  The  gradually  increasing  SIR  for 
positive  distortion  results  from  the  decreasing  autocorrelation  signal  shown  in  plot  (a),  while 
the  crosscorrelation  remains  constant.  Encoded  spectral  patterns  for  which  the  outermost 
features  at  long  and  short  wavelengths  are  complementary  do  not  cause  an  imbalance  at  the 
receiver  for  positive  distortion. 

While  it  might  have  been  expected  that  limited  resolution  would  be  the  critical  parame¬ 
ter,  the  reduction  in  contrast  between  adjacent  spectral  chips  does  not  degrade  the  balance 
of  the  photodiodes.  Equal  amounts  of  light  are  directed  to  the  wrong  detector,  and  the 
crosscorrelation  signal  remains  zero.  However,  the  limited  resolution  does  degrade  the  mag¬ 
nitude  of  the  autocorrelation  signal,  since  light  diffracted  into  neighboring  spectral  chips  will 
be  received  by  the  wrong  photodiode  and  subtracted  from  the  desired  signal  of  the  other 
photodiode.  This  effect  is  apparent  in  the  measurements  of  Fig.  3,  where  the  autocorrela¬ 
tion  peaks  for  code  Jb  are  roughly  half  those  for  codes  Ja  and  Jc-  Although  each  of  the 
super-codes  is  length  32,  the  code  sequence  for  Jb  is  composed  primarily  of  alternating  I’s 
and  O’s,  while  the  sequences  for  Ja  and  Jc  are  composed  primarily  of  repeated  segments  (11 
and  00),  reducing  the  effects  of  limited  resolution. 

To  better  understand  the  origin  of  the  wavelength  scale  distortion  of  the  decoder,  an 
analytical  sensitivity  analysis  was  performed  using  the  plate  factor  of  the  decoder.  The 
plate  factor  P  expresses  the  number  of  nanometers  of  optical  bandwidth  per  millimeter  of 
spatial  extent  in  the  mask  plane.  For  a  grating  with  a  groove  spacing  of  d,  a  collimation 
mirror  of  focal  length  /,  and  an  observation  wavelength  A  at  a  dispersion  angle  of  /5, 

^  dcosp 

A  simple  derivative  of  this  expression  reveals  the  one  to  one  dependence  5P/P  =  —Sf/f. 
Given  the  precision  of  our  breadboard  design  and  the  commercial  tolerances  of  the  mirrors, 
we  are  confident  that  the  distortion  is  not  due  to  a  5%  focusing  error.  Furthermore,  in  the 
process  of  searching  for  the  cause  of  the  wavelength  scale  distortion,  we  observed  that  spectral 
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components  in  the  decoder  were  dispersed  at  steeper  angles  relative  to  the  encoders.  The 
angle  of  incidence  to  the  grating  was  well  controlled  by  the  breadboard  design  at  a  =  63.0 
degrees,  and  the  observation  wavelength  was  precisely  controlled  by  a  spatial  mask  placed 
at  the  encoder  to  create  a  sharp  spectral  feature  at  approximately  1532  nm.  We  measured  a 
shift  in  dispersion  angle  of  5(5  =  +0.76  degrees  from  the  expected  angle  of  =  71.46  degrees, 
ov  5P/P  —  1.06%.  Taking  a  derivative  of  the  plate  factor  P  with  respect  to  the  dispersion 
angle  /?  and  using  the  grating  equation,  we  find 


SP 

P 


'1  —  ((A/d)  —  sina)^  +  tan/? )  d/?. 


(2) 


Evaluating  this  expression  with  the  values  already  mentioned  and  d  =  833.33  nm  gives 
5P/P  =  —0.0422,  which  compares  favorably  with  the  observed  value  of  a  5%  distortion. 
The  minus  sign  from  the  calculation  is  consistent  with  the  fact  that  the  measured  spectrum 
shown  in  Fig.  5  occupies  a  wider  spatial  extent  than  desired,  corresponding  to  a  reduced  plate 
factor  P.  The  certainty  in  the  value  of  the  incident  angle  a  and  the  observation  wavelength 
leads  us  to  examine  the  sensitivity  of  P  to  the  groove  spacing  of  the  diffraction  grating. 
A  5%  distortion  requires  only  5d/d  =  —0.0023!  Thus,  the  observed  distortion  could  result 
from  a  groove  spacing  of  831.41  nm  instead  of  833.33  nm  or  a  groove  density  of  1202.7/  mm 
versus  1200.0/  mm.  While  the  three  holographic  diffraction  gratings  used  in  the  testbed  were 
fabricated  in  the  same  batch,  they  were  written  independently  and  such  small  variations  seem 
possible. 


5  Improved  Correlation  Measurements 

Having  determined  that  the  wavelength  scale  distortion  can  be  a  critical  parameter  in  the 
correlation  process,  we  conducted  an  experiment  to  correct  the  problem.  While  it  is  possible 
to  increase  the  plate  factor  by  increasing  the  incident  angle  a,  which  decreases  13,  this  has 
the  unwanted  effect  of  decreasing  the  return  angle  of  light  from  Mirror  1  while  increasing 
the  return  angle  from  Mirror  2  (see  Fig.  1).  A  simple  alternative  is  to  use  a  mask  pattern 
at  the  decoder  that  compensates  for  the  distortion.  The  pattern  for  super-code  Jc  was 
designed  and  fabricated  with  varying  amounts  of  distortion  (-2.5%,  0.0%,  +2.5%,  +5.0%, 
and  +7.5%)  to  test  the  compensation.  Shown  in  Fig.  8  is  the  autocorrelation  signal  meaured 
with  the  original  mask  at  the  encoder  and  the  intentionally  distorted  mask  at  the  decoder. 
The  received  electrical  signal  is  clearly  maximized  for  the  pattern  with  +5%  distortion  that 
matches  the  expansion  of  the  spectral  pattern.  The  measured  autocorrelation  power  drops 
29.4%  going  from  a  +5%  to  a  zero  distortion  mask.  The  simulated  autocorrelation  curve  for 
code  Jc  resembles  Fig.  7  (a),  with  a  more  gradual  reduction  for  increasing  distortion.  Because 
the  code  sequence  for  Jc  has  a  number  of  repeated  segments,  a  smaller  proportion  of  light 
is  directed  to  the  wrong  detector.  The  simulation  predicted  a  47.6%  drop  in  autocorrelation 
power. 

Measurements  of  improved  crosscorrelation  performance  were  less  conclusive  due  to  am¬ 
biguity  in  the  horizontal  position  of  the  mask,  as  already  discussed.  For  the  decoder  mask 
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Figure  8:  The  autocorrelation  signal  measured  on  the  spectrum  analyzer  as  a  function  of 
the  distortion  of  the  mask  pattern  at  the  decoder.  The  data  clearly  shows  that  the  signal  is 
maximum  at  a  +5%  distortion,  in  agreement  with  the  measured  expansion  of  the  received 
spectral  pattern. 

with  distorted  patterns,  this  problem  was  more  pronounced  because  the  spatial  width  of 
the  patterns  varies.  As  a  result,  it  was  possible  to  make  adjustments  during  the  measure¬ 
ments  such  that  the  crosscorrelation  signal  was  on  average  35  dB  below  the  autocorrelation, 
regardless  of  the  distortion  of  the  decoder  mask.  Nevertheless,  having  demonstrated  the 
improvement  in  autocorrelation  performance  that  our  simulation  results  suggest,  we  believe 
that  similar  improvements  in  crosscorrelation  performance  are  also  being  acheived  for  code 
combinations  with  sensitivity  to  distortion. 


6  Conclusion 

Passive,  all-optical  correlators  are  attractive  for  high-speed  network  applications  where  asyn¬ 
chronous  operation  and  independence  of  data  rate  are  desired.  They  are  able  to  process 
the  extremely  high  combined  data  rate  on  a  common  network  channel  while  the  detection 
electronics  operate  at  only  individual  user  rates.  However,  the  degradation  of  correlator 
performance  has  a  direct  impact  on  the  capacity  of  a  network  to  support  simultaneously 
active  users.  Our  numerical  simulations,  analytical  results,  and  experimental  measurements 
have  identified  wavelength  scale  distortion  as  a  critical  optical  parameter  in  the  correlation 
process  and  illustrated  the  importance  of  well-matched  physical  characteristics  when  using 
multiple  gratings  or  other  despersive  elements.  Correlator  designs  that  allow  dynamic  tuning, 
such  as  arrayed-waveguide  devices  featuring  adjustable  optical  paths  through  temperature 
or  mechanical  tuning,  might  overcome  this  peformance  limitation.  However,  as  the  effects  of 
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distortion  are  reduced,  the  impact  of  limited-resolution  and  other  device  specific  efifects  will 
set  a  new,  but  lower  limitation  on  correlator  performance. 
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Imaging  characteristics  of  poly(methyi  methacrylate)  at  vacuum  ultraviolet 
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We  report  the  depth  versus  exposure  characteristics  of  poly(methyl  methacrylate)  (PMMA)  in  the 
vacuum  ultraviolet  (VUV)  region  from  59  to  128  nm.  Calculated  absorption  coefficients  of  the 
PMMA  are  also  presented.  The  depth  of  the  features  on  the  PMMA  were  in  the  range  of  20-80  nm 
created  by  1-50  mJcm”^  exposures  and  were  measured  with  an  atomic  force  microscope  (AFM). 

We  found  that  the  AFM  can  reliably  measure  depth  differences  of  only  2  nm.  The  sensitivity  of  the 
AFM  permits  exposures  to  be  reduced  for  high  resolution  microscopy  and  holography  even  in  the 
case  of  highly  absorbing  films.  ©  1997  American  Vacuum  Society.  \SUl?>A-l\\yi{91)\2A0A-0\ 


L  INTRODUCTION 

Microscopy  using  contact  printing  or  holography  requires 
a  well  characterized,  high  resolution  recording  film.  One  of 
the  earliest  and  best  characterized  electron  and  x-ray  resists 
is  poly  (methyl  methacrylate)  (PMMA).^“^  In  contrast,  how¬ 
ever,  to  lithography  applications  —  where  the  full  removal 
of  the  resist  is  desirable  in  the  exposed  areas  —  contact 
imaging  and  holography  require  a  resist  with  vertical  (depth) 
resolution.  The  attainability  of  a  “grey-scale”  pattern  in  the 
resist  is  important  in  these  applications  for  reliable  and  re¬ 
producible  mapping  of  the  light  intensity  in  a  wide  exposure 
range.  In  the  soft  x-ray  region,  microscopy  using  PMMA  as 
a  resist  material  is  well  developed.^  In  PMMA,  a  positive 
resist,  bonds  are  broken  by  incident  short  wavelength  radia¬ 
tion,  decreasing  the  average  molecular  weight  of  the  long- 
chain  molecules.  This  increases  the  solubility  of  the  resist  in 
a  suitable  solvent,  resulting  in  a  pattern  of  surface  heights 
corresponding  to  the  exposure.  The  surface  pattern  can  be 
measured  by  a  scanning^  or  transmission  electron 
microscope, but  an  atomic  force  microscope  (AFM)  pro¬ 
vides  more  reliable  quantitative  data.^’^ 

Vacuum  ultraviolet  (VUV)  microscopy  has  the  potential 
to  supplement  x-ray  and  electron  microscopy,  because  the 
contrast  mechanism  is  quite  different.  VUV  radiation  inter¬ 
acts  with  the  outer  valance  electrons  and  molecular  absorp¬ 
tion  is  sensitive  to  composition  and  structure.  In  addition, 
several  reasonably  practical  incoherent  and  coherent  sources 
exist  in  the  80-130  nm  region.  However,  almost  no  studies 
have  been  made  of  PMMA,  or  other  films,  for  VUV  imaging. 
While  there  are  lithography  studies  in  the  110-260  nm 
region, this  information  is  only  partly  relevant  since  the 
requirements  of  lithography  are  quite  different  from  imaging. 
For  lithography,  the  critical  parameter  is  the  exposure  re¬ 
quired  to  completely  remove  the  film;  imaging,  however, 
requires  quantitative  knowledge  of  the  response  function  and 


“^Present  address:  Department  of  Optics  and  Quantum  Electronics,  JATE 
University,  H-6720  Szeged,  Dom  ter  9,  Hungary. 

^^Permanent  address:  Research  Institute  for  Solid  State  Physics,  P.O.  Box 
49,  H-1525,  Budapest,  Hungary. 

‘^^Electronic  mail:  young@rice.edu 


a  reasonable  dynamic  range,  or  intensity  level  resolution. 
The  relationship  between  exposure  and  surface  height  is  gen¬ 
erally  poorly  known  even  for  common  photoresists.  Basic 
and  important  optical  constants,  such  as  absorption  coeffi¬ 
cient  at  VUV  wavelengths  are  also  lacking. 

In  this  work  we  will  present  the  depth-exposure  charac¬ 
teristics  of  PMMA  in  the  59-128  nm  spectral  range,  for 
exposure  values  of  1-50  mJ  cm”^.  The  surface  depths  were 
measured  with  an  AFM  and  were  in  the  range  of  2-80  nm. 
Since  the  AFM  can  reliably  measure  surface  height  differ¬ 
ences  of  only  2  nm,  very  low  exposures  can  be  used  even  at 
wavelengths  that  are  highly  absorbing.  We  have  used  our 
depth  versus  exposure  data  to  calculate  the  VUV  absorption 
coefficients  of  PMMA  using  a  simple  model  for  exposure 
and  dissolution. 

II.  EXPERIMENTAL  DESCRIPTION 

A.  Preparation  and  development  of  the  PMMA  resist 

The  samples  for  exposure  were  12.5  mm  diameter  pol¬ 
ished  glass  disks  spin  coated  with  PMMA.  The  resist  was  a 
commercial  solution  of  4%  PMMA  (molecular  weight  of 
950  000)  in  chlorobenzene.^^  The  top  surface  of  the  substrate 
was  flooded  with  the  PMMA  solution  and  then  spun  at  5000 
rpm  for  25  s.  A  second  layer  of  PMMA  was  spun  over  the 
first  to  provide  a  thicker  film.  The  resulting  film  was 
290±10  nm  thick  as  measured  by  the  AFM.  The  substrate 
and  film  was  then  baked  at  140  °C  for  at  least  30  min  in  a 
vacuum  oven  at  ^^500  mbar  pressure. 

Following  exposure  the  film  was  developed  in  pure  me¬ 
thyl  isobutyl  ketone  for  60  s  at  25  ±0.5  ®C,  during  which 
time  the  developer  was  agitated  mildly.  Following  develop¬ 
ment,  the  sample  was  immediately  rinsed  in  isopropyl  alco¬ 
hol  for  5  s,  dried  with  a  jet  of  clean  freon,  and  baked  again  at 
140  for  30  min. 

B.  System  setup 

Figure  1  shows  the  experimental  arrangement.  Broadband 
incoherent  VUV  light  was  generated  using  a  laser-produced 
plasma.  Pulses  from  a  2 -switched  NdiYAG  laser  were  fo- 


828  J.  Vac.  Sci.  Technol.  B  15(4),  Jul/Aug  1997  0734-211  X/97/15(4)/828/5/$1 0.00 


©1997  American  Vacuum  Society  828 


830 


Ferincz,  Toth,  and  Young:  Characteristics  of  poly(methyl  methacrylate) 


830 


Fig.  4.  An  atomic  force  micrograph  of  a  contact  printed  pattern  on  the 
surface  of  PMMA. 


exposure  was  obtained  by  simply  multiplying  this  average 
exposure  value  per  laser  pulse  by  the  total  number  of  pulses 
in  an  exposure. 

In  order  to  account  for  possible  long-term  changes  in  the 
efficiency  of  the  laser-produced  plasma  and  the  light  collec¬ 
tion  system,  we  interrupted  each  exposure  to  take  source 
spectra  at  Port  No.  2.  The  number  of  spectra  taken  was  var¬ 
ied  to  achieve  good  average  values.  Three  spectra  were  mea¬ 
sured  for  the  shortest  (3  and  4  min)  exposures  and  seven  for 
the  longest  (60  min)  exposures.  The  exposure  energy  was 
calculated  by  averaging  the  spectra. 


D.  Depth  measurement  by  AFM 

The  transverse  and  vertical  scans  of  the  AFM^^  were  cali¬ 
brated  using  a  calibration  grating  and  height  standards.  A 
typical  atomic  force  micrograph  of  a  contact  printed  test  pat¬ 
tern  on  the  surface  of  the  PMMA  is  shown  in  Fig.  4.  To 
obtain  a  good  estimate  of  the  depth  of  a  given  hole,  four 
'--50X50  jim  (256  by  256  pixels)  scans  were  made  using 
different  scanning  directions  and  the  results  averaged.  Typi¬ 
cally,  the  standard  deviation  of  the  pixel  height  from  the 
average  was  3  nm.  However,  because  of  the  large  number 
of  pixels,  this  deviation  did  not  effect  the  depth  measurement 
directly,  since  the  average  height  of  the  pixels  were  com¬ 
pared.  However,  this  high  frequency  noise  did  introduce  er¬ 
rors  during  the  removal  of  the  apparent  surface  curvature. 
Apparent  surface  curvature  is  an  inherent  byproduct  of  the 
AFM  transverse  scanning  mechanism,  and  it  can  be  well 
approximated  by  a  second  order  polynomial.  The  elimination 
of  the  curvature  was  performed  by  fitting  a  second  order 
polynomial  to  each  data  row  of  the  scans.  The  accuracy  of 
this  fit  is  reduced  by  the  pixel  height  variation  because  of  the 
small  number  of  pixels  in  each  row.  This  was  particularly 
true  where  the  data  row  contained  values  from  both  the  ex¬ 
posed  and  unexposed  parts,  since  the  unexposed  part  was  not 
included  in  the  fitting.  We  estimate  that  the  vertical  resolu¬ 
tion  was  limited  to  2  nm  by  this  effect. 


Fig.  5.  Developed  film  depth  vs  exposure  for  PMMA  at  two  different  wave¬ 
lengths.  The  symbols  with  error  bars  are  the  measured  height  differences 
between  exposed  and  unexposed  regions,  while  the  lines  are  the  results  of  a 
model  fitted  to  the  data. 

III.  RESULTS  AND  DISCUSSION 

Measurement  of  exposure  patterns  such  as  Fig.  4  resulted 
in  characteristic  curves  of  depth  versus  exposure,  S,  at  spe¬ 
cific  wavelengths  spaced  4.25  nm  over  the  range  of  59-128 
nm.  Figure  5  shows  two  such  curves;  the  symbols  represent 
measured  data  while  the  lines  are  the  result  of  the  exposure- 
dissolution  model  discussed  below.  This  set  of  curves  can  be 
usefully  summarized  by  defining  a  specific  exposure  flux, 
Sp,  for  each  wavelength  that  will  produce  a  30  nm  deep 
pattern.  Figure  6  shows  a  plot  of  Sp  over  the  wavelength 
range  measured.  While  the  30  nm  criteria  is  arbitrary,  we 
found  that  patterns  of  this  depth  resulted  in  low  noise,  repro¬ 
ducible  measurements.  Also  shown  in  Fig.  6  is  the  calculated 
derivative  of  pattern  depth  versus  exposure,  G  =  dzldS, 
evaluated  at  Sp .  The  parameter  G  is  a  measure  of  the  expo¬ 
sure  resolution;  a  large  G  indicates  there  are  more  distin¬ 
guishable  intensity  levels  within  a  given  exposure  range. 
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Fig.  6.  Specific  exposure  Sp  required  for  a  developed  depth  of  30  nm, 
(triangles),  as  a  function  of  wavelength,  and  G,  the  derivative  of  depth  with 
respect  to  exposure  at  Sp  (filled  squares),  characterizing  the  possible  grey¬ 
scale  levels  for  imaging  applications. 
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reconstruction  of  contact  and  holographic  recordings.  We 
found  that  exposures  as  snmall  as  2  mJcm“^  can  produce 
patterns  20  nm  deep.  Since  the  AFM  can  resolve  '--2  nm 
depth  differences,  about  ten  vertical  or  exposure  levels  are 
distinguishable.  This  exposure  level  is  significantly  lower 
than  that  previously  used  with  other  depth  measurement 
techniques. Further  improvement  should  be  possible  by  de¬ 
creasing  the  pixel  noise,  perhaps  by  modifying  the  film 
preparation  and  developing  procedures. 
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